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General Conclusion . 

Dr. Kiistner’s magnitude equation as determined by a com¬ 
parison with Oxford photographic measures is in substantial 
agreement both in B,.A. and Decl. with that found by him, 
using the method of screens ; but the quantities are too small to 
settle the question whether the screen method is apt to give 
values which should be multiplied by a factor, as suggested in a 
former paper. 


On Periodic Orbits in the Neighbourhood of Centres of Libration. 

By H. C. Plummer, M.A. 

i. It is well known that there are five exact solutions of the 
problem of three bodies in each of which the bodies preserve an 
unvarying configuration which revolves with a uniform velocity. 
It is also known that when the third body is of infinitesimal 
mass compared with the other two, it can describe small periodic 
orbits in the vicinity of the points where exact solutions exist. 
These points have been given by Gylden the name of centres of 
libration , and a body describing an orbit of the kind indicated 
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Nov. 1901. Neighbourhood of Centres of Libration. 7 

has been called by Professor Darwin an oscillating satellite . 
Such orbits have been investigated recently by Dr. Charlier in a 
very interesting and suggestive paper.* The method which he 
has adopted is extremely simple and familiar, and several of his 
results had been given previously by Professor G. H. Darwin, in 
whose work the subject only presented itself incidentally, and by 
Mr. Moulton in his paper on the origin of the Gegenschein.\ I 
wish here to discuss some of Dr. Charlier’s results in a slightly 
more general manner. It is true that the direction in which 
greater generality can be secured without increased difficulty 
will lead us away from the problem of three bodies. But it is 
possible that the conclusions which follow will prove interesting 
in themselves, that they may conceivably find an application to 
some physical problem of a totally different order, and that 
something more may be learnt with regard to their general 
character and limitations. At any rate any attempt to draw 
wider attention to a subject which presents such a promising 
field for investigation may perhaps be justified. 

2. Let us consider the motion of a particle in a rotating field 
of force. The forces are to be derived from a complete potential 
function O, the precise form of which need not at present be 
defined except as some uniform function of x and y , the 
coordinates referred to rotating axes. Then the equations of 
motion of the particle may by a suitable choice of units be 
expressed in the form 

dO, 

*- 2 y=d* 

3 Qi 

y+2±=^ 


As is well known these apply in particular to the motion of a 
body of infinitesimal mass under the influence of two bodies 
revolving in circles about their centre of gravity. Now let a 
centre of libration be defined as a point whose coordinates (x, y ) 
satisfy the relations 

3a__3o,__ 

dx 3 y 

At such a point equilibrium, which may be either stable or 
unstable, is possible. Now let a surface be defined by the 
equation z=il, and let the ordinary notation be adopted, viz. 

3 z 3 z 

Tx =P 0y = ? 

d 2 z _ d 2 z __ &z_ 

dx 2 dxdy S dy 2 

* Meddelanden fran Lunds Astronomiska Observatorium , No. 18. 
t Astr. Journ ., No. 483. 
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Now if in the equations of motion we substitute cr-j-£, y-\-rj for 
x and y and consider ( x , y) a fixed point, we obtain 

£ . 30 . f 3 2 12 . 9 2 & 

2ri dx dx 2 ** dxdy 

, L dn , y d 2 £l , d 2 Q, 

n + 2 i = T y + ‘ + v dr + " 

Hence if (x, y) is a centre of libration so that p=q-=zo and if we 
neglect second and higher powers of l and rj we may write 

l— 2fi=r£ + 8r) 
rj + 2£=s£ + £7 

Then £ =Ae xt , rj=Be Xi satisfy the equations provided 

A(\ 2 — r) —B(2 X + s) = o 
A( 2 \—s) -f- B(X 2 —£)=o 

i.e. (\ 2 —r)(\ 2 —t) + 4\ 2 — s 2 =o 

or X 4 + A. 2 (4—r — t) -\-rt — s 2 =o 

Stability depends on the nature of the roots of this equation. 

3. Now let p 19 p 2 be the principal radii of curvature at the 
corresponding point of the surface z = O. Since p = q = o they 
are the roots of the quadratic 

{rt—s 2 )p 2 + (r + t)p+ 1 =0 

Hence we may rewrite the preceding equation in X 2 in the form 


X 4 + A 44 + ! + -)+— = 
X Pi P2 y Pl°- 


P*' 

By means of this equation the conditions of stable equilibrium 
at a centre of libration are immediately connected with the pro¬ 
perties of the surface z = O. For perfect stability it is necessary 
that both values of X 2 shall be real and negative. Hence we 
must in the first place have both p x and p 2 of the same sign, the 
surface synclastic at the point, and z a real maximum or a real 
minimum. The second condition is that 

4 +- + —> - - 

Pi P 2 + >/ p I j0 2 

If z is a maximum, p x and p 2 are positive, and the condition 
is necessarily fulfilled. If z is a minimum, p 2 and p 2 are negative, 
and the condition may be written 


-= + 


+ \/—pi + V — P2 


.<2 
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Nov. 1901. Neighbourhood of Centres of Libration. 9 

Thus, while in a stationary field of force equilibrium is stable 
where the force function is a maximum only, in a rotating field 
the function £2, which it must be remembered includes the rota¬ 
tion potential, may be either a maximum or a minimum, provided 
in the latter case a certain inequality is satisfied. On the other 
hand we may have one value of X 2 negative and the other posi¬ 
tive. This corresponds to curvatures of opposite sign or an 
anticlastic point on the 12 surface. In this case a mode of 
periodic motion is combined with an unstable mode which must 
sooner or later prove fatal to the stable form. Neither purely 
positive nor imaginary values of X 2 can be admitted in a search for 
periodic orbits. 

4. I consider first the far more interesting case in which two 
periodic solutions are possible. Let us assume 

£=A cos (mt—a) v=k cos (mt—f3) 

as a periodic solution of the equations 

1 —2 rj=r£ — sr) 

The conditions obtained by substitution are 

(r+m 2 )h cos a+sk cos /?= — 2mk sin ft 
(r + m 2 ) h sin a + sk sin 2mk cos /3 

(t+m 2 ) k cos /3 + sh cos a= 2mA sin a 
(t -f m 2 ) k sin /? + sh sin a = — 2 mh cos a 

Hence we obtain directly 

(r + m 2 ) h sin (a — 

, sk sin (a — j3) = — 2 mk cos (a—ft) 

(t + m 2 ) A sin (a— (3)— 2 mh 

sh sin (a— fi)= — 2 mh cos (a—/ 3 ) 

which are equivalent to 

A 2 (r + m 2 )= k 2 (t + m 2 ) 
s tan (a—/ 3 )= — 2m 

4m 2 cosec 2 (a— /3) = (r + m 2 ) (t + m 2 ) =4m 2 + s 2 

Hence the periodic orbit is obtained by combining two simple 
harmonic motions at right angles. These motions are charac¬ 
terised by a constant ratio of the amplitudes and a constant 
difference of phase for each centre of libration and assigned 
period. The orbits may be illustrated by Lissajou’s figures for 
tuning-forks in unison. This is a very simple way of stating the 
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connection between the integration constants. The true arbi¬ 
trages are, of course, the epoch and the scale of the orbit. 
Except for these the character of the orbits is absolutely deter¬ 
mined by the nature of the field of force. 

5. The equation of the orbit can be written down at once in 
the form 

f + S-§ cos(a “^ )= ^ 2(a_/3) 

It will be found that an important distinction must be drawn 
in what follows, according as equilibrium corresponds to a 
maximum or minimum value of Ll . The latter case may be con¬ 
sidered first. Then both r and t are positive, and we may put 

c 2 =h 2 (r + m 2 )=k 2 (t + m 2 ) 

and since 

s(r-\-m 2 )hk sin 2 (a—/?)= — ^m 2 k 2 cos (a—/ 3 ) 
we have also 

— hie . s sec (a— p)—c 2 

The orbit becomes 

(r + m 2 )£ 2 + (t + m 2 )r \ 2 + 2 =4m 2 c 2 / (4m 2 + s 2 ) 

and the axes and foci are given by 

5 2 —if hf c 2 

r—t s 4m 2 -fs 2 

The directions of the axes depend, then, on the form of the 
potential alone, and the major axis lies in the first and third 
quadrants if s is negative, in the second and fourth if s is 
positive. But for the surface z = O 

3 + fo= 0 +p£ + qrj +i(r£ 2 + 2 s£rj + trf) +. . 

or 2(k;=r£ 2 + 2 S%r) + tr ) 2 

is the indicatrix corresponding to a centre of libration. The 
axes and foci are given by 

£ 2 — 7] 2 _£77_(2 lz ) 2 

r—t s s 2 —rt 

and s 2 —rt=—p I ~ 1 p 2 ~ 1 - is negative. Hence the major axes of the 
orbit and of the indicatrix are parallel. An alternative form of 
statement is that the axes of the orbit are parallel to the 
tangents to the lines of curvature at the corresponding point 
of the O surface. 
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6. The semi-axes a , b of the orbit are given by 

- + * =(r+t+2m 2 ) (4m 2 + s 2 ) /4 m 2 c 2 
Oj 2 o 2 

— 4 -=(4m 2 + s 2 ) 2 / 4^ 2 c 4 
a 2 & 2 v 7 


Hence a 2 and 6 2 are the roots of the quadratic 

(4m 2 + s 2 ) 2 x 2 —(r+t+2771 2 ) (4m 2 + s 2 )c 2 a? + 4m 2 c4=o 
2 (4m 2 + s 2 ) . x/c 2 =r + t + 2 m 2 

+ s/ {4m4 + 4m 2 (r + 1 —4) + (r + t) 2 } 

=r + £ -f- 2 to 2 +n/ {4(s 2 — r£) + (r + £) 2 } 

= 2 m 2 - ( - + ± f- --) 

\Pi P2J \P 1 P2/ 

We may therefore write 

a^—c 9 [m 9 ——^ /(4m 2 + s 2 ) 

b 2 =c 2 ^m 2 -/ (4m 2 + s 2 ) 


where p.,. is the numerically smaller of the two radii of curvature, 
for these are both negative. Accordingly the eccentricity e is 
given by 

e 2 =(p 2 -'-p-')/(m 2 -pr T ) 

and we notice that a circular orbit corresponds to an umbilic on 
the O surface. 

7. We can now examine the initial conditions and the relation 
of our constant c to the constant C of the Jacobian integral 

W= 2 0-C 


Let (£ 0 . rj 0 ) be the position at the time t= o. Then since 
£=h cos (mt —a) i]—k cos (mt — ( 3 ) 

we have 

? 0 =h cos a Vo—k cos f 3 

£ 0 =mh sin a rj 0 —mk sin (3 
.\ Y 2 =? 0 2 +7) 0 2 =m 2 (A 2 + k 2 ) - ra 2 (£ Q 2 + rjo 2 ) 

But 

20 , = 2O 0 + 2 + trio 2 + 2 s£ 0 r ] 0 

where £l Q is the value of £1 corresponding to the centre oi 
libration. If 2 0 0 =C 0 

C=C 0 -m 2 (h 2 +k 2 ) + (r + m 2 )& + (t+m 2 )r ) 0 2 + 2 s| 0 r/ 0 
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Thus for a given value of C—C OJ £ 0 and rj Q must satisfy the 
equation of a certain ellipse which can be no other than the orbit 
already found. Hence by comparison we find 

C—C 0 =4m 2 c 2 /(4m 2 + s 2 )—m 2 (h 2 4- k 2 ) 

_4m 2 c 2 _m 2 c 2 (2m 2 +r -M) 

4m 2 + s 2 (m 2 4- r)(m 2 + 1 ) 

=m 2 c 2 (4 — r—t— 2m 2 ) / (4m 2 + s 2 ) 

Now if m r 2 and m 2 2 are the roots of the quadratic giving m 2 


or 


m I 2 +m 2 2 =4— r—t 
(2 m 2 +r 4- 1 — 4) + (2m, 2 + r + t— 4)=o 


Hence C—C Q clearly has opposite signs according to the 
period selected, that is, one period and one eccentricity corre¬ 
spond to an increase, and the other period and eccentricity to a 
decrease of the energy constant as compared with the value 
required for relative rest. Thus when stable equilibrium exists 
at a point where 11 is a minimum a given value of C determines 
the orbit uniquely. Dr. Charlier has noticed this in a particular 
case, but his proof is not of a general character. Another 
peculiarity is that one family of orbits is bounded by a set of 
“ limiting curves,” while the other is not. 

8. The period of the orbit is 2 tt /ra. But the units have been 
so chosen that the period of rotation of the field of force is 27 r. 
Therefore m may be defined as the ratio of the period of the 
rotating field to the synodic period of the orbit. The direction 
of motion is determined by the sign of 

£rj——mhk sin (a— j 3 ) 

= — 2 m 2 c 2 / (4m 2 + s 2 ) — —mob 


The angular momentum being negative, the orbits of both families 
are retrograde. If it is desired to write down the equation of an 
orbit in a form involving only the energy constant and quantities 
depending only on H, we are led to 

(r + m 2 )f 2 + (t + m 2 )?? 2 + 2^77=4(C — 0 o ) / (4 — r—t —2m 2 ) 

=4%(C — C 0 ) 


where 


n~ 2 =(4 ^— r — t ) 2 — 4(r t—s 2 ) 

= W-4--) 4-i6 

\Pi Pa/ \Pi Pa/ 


It may be observed that it is the family of shorter period and 
smaller eccentricity which requires C < C Q . 

9. Throughout §§ 5-8 it has been assumed that 12 is an 
absolute minimum at the centre of libration, and that accordingly 
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r and t are positive. This indeed happens to be the more inter¬ 
esting case in so far as the problem of three bodies is concerned 
because there the function 12 does not as a matter of fact possess 
a truly maximum value. But if 12 is an absolute maximum 
both r and t are negative. It will be convenient to denote the 
quadratic giving m 2 by y(m 2 )=o. By substitution we obtain : 


/ (+ 00 ) is positive. 


are negative. 


/(—r)=4r—s 2 
f(o) =rt — s 2 is positive. 


. Hence — r and — t separate the two positive roots m 2 and 
m 2 2 . It is clear then that when li is a maximum mf + r and 
m 2 + t are negative, while m 2 2 +r and m 2 2 + t are positive if 
m 2 2 >m I 2 . Further, since 


4— r — t — 2m T 2 > o 

we must have 

4 —r — t — 2m 2 < o 

The orbits of shorter period, 2n are then similar to the orbits 
investigated in the case when 12 is a minimum, and they require 
that C < C c . It is unnecessary to consider this family in fur¬ 
ther detail. It is only in the expressions for the semi-axes and 
eccentricity given in § 6 that it is necessary to interchange the. 
two radii of curvature on account of the change of sign. 

10. The family of orbits characterised by the longer period 
27 r/m x possesses features, however, which require separate 
notice. In the place of +c 2 in the foregoing it is necessary to 
write — c 2 . The formula given for the foci in § 5 shows that the 
major axis of a member of this family is at right angles to the 
corresponding axis of a short-period orbit. Further from § 8 it 
appears that the direction of motion is reversed and orbits of 
this family are direct. Since both c 2 and (4 — r — t — 2m 2 ) have 
opposite signs in the two families we again have C < C D . Con¬ 
sequently to any given value of C C G one orbit of each family 
corresponds. Indeed the Jacobian integral suffices to prove that 
C j> C Q ; for otherwise since 2l2 Q ^212, the velocity of the displaced 
body would be imaginary. Hence the extension of Dr. Charlier’s 
statement given in § 7 of this paper, must be restricted to cases 
where 12 is a minimum. The expressions given for the semi¬ 
axes in § 6 require no alteration, although the radii of curvature 
are positive : for (m x 2 —p x -1 ) and (ra x 2 —■ p 2 -1 ) must be negative, 
and therefore p z ~ x > p 2 " J . The most interesting features of the 
family of orbits considered here are undoubtedly the change from 
retrograde to direct motion and the rotation of the major axis 
through a right angle. 

n. When 12 is a maximum at the centre of libration, we 
have seen that two orbits correspond to one value of C. Again, 
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through a point two orbits can be drawn, and for these the 
values of C will in general be different. There are, however, 
points where two orbits of different families, and the same C 
meet. The orbits will be 

(m x 2 + r )? 2 + (m/ + t)rj 2 + 2S?r/= + 4™(C—C Q ) 

(ra 2 2 +r )? 2 + (^ 2 2 + t)n 2 + 25???= — 4^(C — C 0 ) 

Hence the locus of their intersection is 

( m ^+^ 2 2 + 2 r )? 2 + ( m x 2 + m 2 2 + 2^)?? 2 + 4 s ? t?=o 

Or since m I 2 + m 2 2 =4 — r — t 

(4 + r—2)£ 2 + ( 4 -r + £)*? 2 + 4 ^=° 

This represents a pair of straight lines equally inclined to the 
axes of the orbit, and this is clear from the geometrical 
interpretation. For since C is the same for both orbits, the 
velocity in each at the point of intersection is the same. Hence 
the central perpendiculars on the tangents are equal, and it 
becomes evident that the above equation represents the locus of 
points where orbits of the two families intersect each other, 
and are equally inclined to the common radius vector. The 
locus will not be real unless 

4s 2 + (r — £) 2 > 16 

i.e . |Oi -I ^p 2 _I > 4 

And unless this condition is satisfied corresponding members 
of the two families will not intersect, so that one orbit will lie 
wholly within the other. 

12. When at a centre of libration the function O is neither 
an absolute maximum nor an absolute minimum, oscillating 
satellites may still exist. In this case only one class of orbit is 
possible, and even this cannot be permanently stable. The 
equation f[m 2 )=-o has one positive root and one negative. Since 
the equation 

( 2m z 2 + r + 1 —4) + (2m 2 2 + r + 1 —4)=o 

still holds, the first expression on the left hand must be positive if 
m/ > o > m 2 2 . Hence (2m r 2 + r + 1 ) is positive, and since (m x 2 + r) 
and (mf + 1 ) must be of the same sign, if an elliptic orbit exists 
they must both be positive. Hence C <r C Q , and the orbits possess 
the same general features as those of shorter period when 12 is a 
minimum. An unstable form of motion arises immediately if 
C > C OJ and this would be the case if motion in the periodic orbit 
be unduly retarded. We should infer that excess rather than 
defect of relative velocity would be favourable to the duration of 
quasi-stable motion, and this may have some bearing on the 
theory of the detention of meteorites in the neighbourhood of a 
centre of libration. We have seen that in general the direction 
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of the major axis of an orbit depends on the sign of s. An 
exception, of course, arises when s=o. This point happens to 
occur in the application of this paragraph, and may properly be 
considered in this connection. The formula given in § 5 for the 
foci of the orbit shows that the axes are parallel to the rotating 
axes of coordinates, and that the major axis is parallel or perpen¬ 
dicular to the line joining the controlling bodies, according as 
r—t< or >0. The only exception to this rule is in the case of 
long-period orbits corresponding to 12 a maximum, when the pre¬ 
ceding inequalities are reversed. 

13. About the special problem which suggested the foregoing 
slightly more general investigation, little need be said here, be¬ 
cause it has been treated with some fullness by Dr. Charlier. 
But as I have had occasion to allude to the properties of the sur¬ 
face defined by the equation 3 = 12, it may not be out of place to 
consider briefly the nature of the surface determined by two 
bodies revolving in circles about their centre of gravity. The 
sum of the masses fi and v is taken as the unit of mass, while the 
unit of time is such that the gravitational constant is unity and 
the unit of length the distance between the masses. Then the 
period T= 27 r, the angular velocity is unity, and we have 

2O^(r I -+j0+v(r,-+jL) 

where r„ r 2 are the distances from jjl and v respectively. Now 
Professor Darwin in his memoir * has drawn a typical series of 
the curves 2 12=const., and if these be regarded not merely as 
a family of curves with an arbitrary parameter, but as a set of 
contour lines of a surface, we easily get a fair notion of the 
general character of that surface. Much indeed can be gathered 
from the simplest reasoning. At the two points where r x =r 2 — 1, 
it is easy to prove that 2 12 is a true minimum and equal to 3. 
The familiar topographical illustration will be useful here. At 
each of these minimum points there is a valley. The function 12 
becomes infinite at the points where the masses are placed, but for 
the present purpose these points may be represented by the 
summits of mountains. The contour lines, which are evidently 
symmetrical with respect to the line joining the masses, must 
rapidly approach a circular form as they recede from the centres 
of force and indicate rising ground in all directions outwards. 
The additional features of the surface must clearly be three 
passes, one running between the mountains and connecting the 
valleys, and two others beyond the two mountains and separating 
them from the ground rising in the distance. Each of the three 
passes has its highest point on the line of symmetry, and these 
three points correspond to the three collinear centres of libration. 
Here the 12 surface is anticlastic, and, as indeed is evident on 
dynamical grounds, stable equilibrium is out of the question. 

* Acta Math ., t. xxi. p. 113. 
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14. Although there cannot be any doubt as to the nature of 
these points, yet it may be permissible to consider them from 
another point of view. In what follows the origin is taken at the 
centre of gravity of the revolving bodies, the latter always lying 
on the rotating axis of x. Since 


we have 


3o 

dr 


2 a = +~}+y (rs + y 

1 ( I\ 9 Q / 


h) 


or 


Also for the collinear centres of libration 

30 

y =°> 35 = ° 

30 dr 1 30 3 r 2 _ 

3 r x ’ dx dr 2 dx 

Moreover 30 _ 30 dr T . 30 3 r 2 

3 y dr x ’ 3 y dr 2 ’ dy 

_ y 0O_j y_ 30 

r x * 3 r x r 2 * 0r 2 

Hence when y = o 


^ 3 2 0 1 30 1 30 

dy 2 r x * 3 r x r 2 dr 2 

0 2 O 

We have now to treat separately the cases which arise. 
(1) Between the masses y and v we have 

dr dr 

r I +r,= i;-p+p=o 

OX OX 

. da_da 

r)r x or . 2 

Hence clearly r x > r 2 if y > v. In this case 

and is negative, for r x < 1. 
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(2) Beyond the mass p we have 

# 0r T- _0r 2 

^2-^1 + I j T: 3 

da? ox 

. an L an 

, . -— 4- —=0 

a r 1 dr 2 

••• 


17 


Hence clearly r z < 1. 


In this case 



and is negative because r T < 1 <r 2 , and therefore the first factor 
is positive and the last negative. The proof that t is negative 
in the case of the centre of libration lying beyond the mas 3 v can 
be inferred from the latter, as there is no restriction on the ratio 
of the masses fi and v . How we have 
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which is necessarily positive. Hence always for the three points 
considered 

rt—s*<o ' 


and therefore the two roots of the quadratic in m 2 are of opposite 
signs. I have thought it worth while to give a definite proof 
of this fact because on this point Dr. Charlier says : “ That is 
what I have found through numerical calculation for different 
values of \i , an algebraic demonstration of this fact seeming to 
be somewhat complicated/' 
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A New Method of Interpolation. By T. C. Hudson, B.A. 
1. Bessel's interpolation formula is 
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2. As it stands this formula is partially symmetrical with 
respect to the middle of the interval ; but there is no symmetry 
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